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ABSTRACT 



We describe singularities in the distribution of polarized intensity as a function of Faraday depth (i.e. the Faraday spectrum) caused 
by line-of-sight (LOS) magnetic field reversals. We call these features Faraday caustics because of their similarity to optical caustics. 
They appear as shaiply peaked and asymmetric profiles in the Faraday spectrum, that have a tail that extends to one side. The 
direction in which the tail extends depends on the way in which the LOS magnetic field reversal occurs (either changing from 
oncoming to retreating or vice versa). We describe how Faraday caustics will form three-dimensional surfaces that relate to boundaries 
between regions where the LOS magnetic field has opposite polarity. We present examples from simulations of the predicted polarized 
synchrotron emission from the Milky Way. We derive either the probability or luminosity distribution of Faraday caustics produced in 
a Gaussian magnetic field distribution as a function of their strength, ¥, and find that for strong Faraday caustics P(f) oc T~ 3 . If fully 
resolved, this distribution is also shown to depend on the Taylor microscale, which relates to the largest scale over which dissipation 
is important in a turbulent flow. 
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1. Introduction 

The availability of broadband receivers in the next generation 
of radio telescopes, such as the Expanded Very Large Array 
(EVLA), the Square Kilometer Array (SKA), and its pathfinder 
arrays such as the Low-Frequency Array (LOFAR), has led 
to the development of a novel technique for imaging polar- 
ized radio emission known a s rotation measure (RM) synthesis 
dBrentiens & de Bruvnll2 005). Derived from a technique origi- 
nally proposed by Burn (I1966I) . RM synthesis allows for the sep- 
aration of polarized sources along the line of sight (LOS) by re- 
constructing the distribution of polarized emission as a function 
of Faraday depth (in the simplest case, Faraday depth is equiva- 
lent to RM). The technique provides several benefits: sensitivity 
is improved by combining measurements from hundreds or even 
thousands of channels over a broad frequency range, the rotation 
measure of point sources can be measured more accurately and 
unambiguously, and the intrinsic polarized emission from differ- 
ent sources along a line of sight can be studied independently. 
Since the Faraday depth depends on the integrated LOS mag- 
netic field, RM synthesis also promises to be an important tool 
for studying magnetism. 

To date, there have been few applications of RM syn- 
thesis, although interest is rapidly increasing as new radio 
telescopes are being commissioned. Its applications have in- 
cluded studies_of_tJie_dffiise polari z ed emission in th e Perseus 



field (Ide Bruvn & Brentiensl l200l iBrenti 
associated with Abell 2255 (iPizzo et al 
the polarized emission in near by galaxies 



l2010h and that 
I2010I) . analysis of 
the WSRT- 



emission in near by galaxies in 
SINGS survey dHealdet al.l l2009h . and the detection of a shell 
of c ompressed magnetic field surrounding a local HI bub- 
ble dWolleben et al.l l201dt) . The RM synthesis technique will 



play a critical role in several upcom ing polari zation surveys, 
e.g. P OSSUM dGaensler et al.ll2010h . GMIMS dWolleben et ail 
I2009I) . and future surveys with LOFAR. 

Th e technique of RM synthesis builds on the work of lBurnl 
(1 1966b who showed that the complex polarized intensity, P, as 
a function of wavelength-squared (A 2 ), is related to the intrinsic 
complex polarized intensity as a function of Faraday depth, F, 
by a Fourier transformation 



P(A 2 ) 



F(cf>, A 2 )e 2i f*d(p. 



(1) 



The polarized intensity is P — Q + iU, where Q and U are the 
Stokes parameters that describe the linear polarization state of 
the emission. The coordinate 0, known as the Faraday depth, is 
a measure of the amount of rotation suffered by a linearly polar- 
ized wavefront as it passes through a magneto-ionic medium. It 
is measured in rad irT 2 and is given by 



4> = a o 



I n, 

Jo 



(z')B 3 (z')dz' 



(2) 



where n e is the number density of thermal electrons, z is distance 
along the LOS, and Bi(z') is the component of the magnetic field 
along the LOS. The constant term is ao = <? 3 l(2mn 2 c A ) = 0.81 
cm 3 (//G pc)~', where e and m e are the charge and mass of the 
electron, respectively. 

The quantity F(<p, A 2 ) is the intrinsic polarized intensity as 
a function of Faraday depth. We refer to this quantity as the 
Faraday spectrum. It is also co mmon l y refer red to as the Faraday 
dispersion function following iBurnl (11966b . Though F generi- 
cally depends on the wavelength, we will assume a separable 
form, F(<p, A 2 ) = f((p)s(A 2 ), allowing for the removal of s(A 2 ) 
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from the integral in Eq. ([T). We will henceforth assume that 
P(A 2 ) can be normalized by the spectral dependence such that 
P(A 2 ) = j Fffle^df. 

The Faraday spectrum is a three dimensional description 
of the polarized emission and may be reconstructed by invert- 
ing Eq. ([TJ. In general, this is not possible because of our in- 
ability to completely sample P for all values of A 2 , but this 
complication is resolved usin g the RM synthesis technique of 
iBrentjens & de Bruvnl (|2005) where the inversion is treated in 
much the same way as it is in synthesis imaging. The incomplete 
sampling of P(A 2 ) is represented by introducing a sampling func- 
tion. The reconstructed Faraday spectrum is then a convolution 
of the true Faraday spectrum and the Fourier transform of the 
sampling function, known as the RM spread function (RMSF). 
The RMSF is a point spread function in (f>, and is analogous to 
the "dirty beam" in synthesis imaging. 

In this paper we are concerned with the properties of the 
Faraday spectrum of extended synchrotron emitting and Faraday 
rotating sources. In particular, we describe singularities in the 
Faraday spectrum that we call Faraday caustics. These arise 
when the LOS magnetic field changes direction, causing a pile- 
up of emission at a single Faraday depth. They are caustics, 
much like the well known optical phenomenon, in the sense that 
they are the result of the non-monotonic mapping of the polar- 
ized emissivity from physical depth space into Faraday depth 
space. Optical caustics, for instance those seen at the bottom 
of a swimming pool on a sunny day, are also caused by a non- 
monotonic mapping of the surface radiation field to its lensed 
image at the bottom. 

Before we derive a mathematical description of Faraday 
caustics, it is helpful to develop a qualitative understanding of 
why they occur. We consider a volume filled with magnetic 
fields, and both relativistic and thermal electrons. This volume 
will produce synchrotron emission and rotate the plane of po- 
larization of a propagating plane wave. The polarized emission 
from each location along a LOS generally maps to a different 
location in Faraday depth space. Therefore we expect to find ex- 
tended structures in RM synthesis observations of such a source. 
When the LOS magnetic field component is changing sign, what 
happens near the point where the field is completely in the plane 
of the sky? In this case, the LOS field is close to zero and the 
Faraday depth will change very little with physical depth be- 
cause, according to Eq. (0, the rate of change in Faraday depth 
with z varies linearly with the LOS magnetic field component. 
This means that all of the polarized emission originating from 
this region will pile up around a single Faraday depth value and 
we expect the resulting Faraday spectrum to be strongly peaked. 

In this paper, we have two aims. First, we wish to in- 
form observers about the fundamental properties of Faraday 
caustics. We expect them to be ubiquitous in the Faraday 
spectra of extended polarized sources, hence they should be 
taken into account when interpreting observational data. We 
first describe the basic features of Faraday caustics using a 
simple description of a magnetic field reversal, and provide 
examples of how they appear in simulations of the Galactic 
synchrotron emission. Second we investigate how observations 
of Faraday caustics can be exploited to uncover information 
about the structure and statistical properties of the underlying 
magnetic field. We present a few ideas and examples. This 
second intention is in line with previous work on the extraction 
of statistical i nforma t ion of magne ti c fields fr om observables 



(12007b. 'letcher & Shukurovl (120071). Waelkens et al ]J2009I) . 
lJunklewitz & Enfilinl ( 12010^ 7 and lOppermann et al .1 (120101) . 

In Sect. |2] we derive the properties of Faraday caustics by 
considering a simple magnetic field near the location where it 
is aligned in the plane of the sky. We examine the observational 
properties of caustics in Sect. [3] presenting examples of simu- 
lated synchrotron emission from the plane of the Milky Way. 
The luminosity function of caustics generated in a Gaussian ran- 
dom magnetic field as a function of their strength is computed in 
Sect. [4] Finally, in Sect. [5] we present our conclusions. 



2. Properties of Faraday caustics 

Faraday caustics appear where the magnetic field lies entirely in 
the plane of the sky, i.e. where B z — 0, in a mixed synchrotron 
emitting and Faraday rotating medium. To characterize the basic 
features of Faraday caustics, we assume the simplest magnetic 
field configuration possible and focus on the region immediately 
surrounding the location where the LOS component is zero. We 
label this point zo, where B?,(zo) = 0, and henceforth for simplic- 
ity we set zo — (the observer is located at z < 0). Ultimately, 
we wish to solve for the Faraday spectrum, F{<p). 

The Faraday spectrum is related to the intrinsic polarized in- 
tensity, P(z), by 



5{<p - <Kz))P(z)dz. 

oo 

The polarized intensity per unit length along z is 

P( z )dz = p(a)I(z)e 2ix(z) dz 

= a s (a,z)\B x (z)\ 1+a e 2ixiz) dz. 



(3) 



(4) 



To simplify the following expressions, we encapsulate all 
non-essential terms into a s . This parameter therefore con- 
tains all of the constant and frequency-dependent terms in 
the expression for the sy nchrotron emissivity, I{z) (see e.g. 
iRvbicki & Lightmanl2004l) . as well as the fractional polarization 
of synchrotron emission, p(a) — (a+ l)/(a + 5/3). It varies with 
Z because it also incorporates the dependence on the cosmic-ray 
electron density. In general, this will be essential for understand- 
ing the detailed structure of a Faraday caustic profile, but in this 
illustrative case we assume that the cosmic -ray electron density 
is roughly constant in the region around zo = 0. The parameter 
a is the spectral index, where /(v) oc v~ a . For convenience, we 
assume that a - I throughout the remainder of this paper. This 
assumption is approximately correct for many optically thin syn- 
chrotron sources and suits our purpose of providing a basic de- 
scription. The angle x is the intrinsic position angle and B ± is the 
magnetic field in the plane of the sky. To compute F, we require 
a description of the z dependence of both B 3 and B ± . 

We restrict our attention to the region immediately surround- 
ing zo = and expand the LOS magnetic field around zo to first 
order with respect to z as 



B 3 (z) = B'(0)z + 



(5) 



where the prime indicates differentiation in z- We recall that 
Z?3(0) = 0, by definition. In this way, we capture the essential 
properties of a field reversal without having to specify a func- 
tional form for B^{z). From Eq. (0, we find that 



by ISpangleri d 19821 1 1983b. lEilekl dl 989albb . lEnfilin & Vogt 
(12003b . iKahniashvili & Vachaspat] (120061) . iMizeva etal 



<p{z) ~4>o = 



aiB'(0)z 2 



(6) 
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where a\ = aon e and we assume that n e (z) is approximately 
constant over this limited range in z- Our illustrative example 
of course becomes inaccurate for distances over which n e starts 
to vary significantly. Such a variation will not affect the basic 
properties of a caustic, so we ignore it for the moment. 

At this point, we can get an idea about the general form of 
F((p). The last expression tells us how the polarized emission 
originating at a given z location maps to 0-space; adopting the 
linear approximation for the z-component of the magnetic field 
given by Eq. [5] we can see that (p is quadratic in z- The value 
(p(z) = <f>o, which occurs when z = Zo = 0, is either a maxi- 
mum or minimum depending on the sign of B' 3 (Q). Therefore, 
the Faraday spectrum will extend only to one side or the other 
of 0o . Each ^-location on this side of <po will be mapped to the 
two physical locations given by the multivalued inverse, z(<p). In 
addition, as z — * the change in Faraday depth, d<p/dz, becomes 
small and the emission over a potentially large range of z values 
"piles up" over a small range of <p values. As a result, we expect 
that F((p) will sharply increase near <po and appear as an asym- 
metric "spike". Smaller values of Is^l lead to a flatter (f>(z ~ 0), 
hence we then expect to find stronger spikes when the LOS mag- 
netic field changes direction over larger distances. 

Since we are primarily concerned with describing the effect 
of a LOS magnetic field reversal, we assume the simplest pos- 
sible form of B ± (z). First, we assume that the orientation of B ± 
is fixed over the limited range of z that we are interested in and 
that it is aligned along the jc-axis. Our magnetic field is then 



B = B ± (x, z)x + B' 3 (0)z z. 



(7) 



To determine B ± (x, z), we use the condition that the magnetic 
field is divergence-free, which leads to 



B ± (x,z) = B ± (0,z)-BU0)x. 



(8) 



With no loss of generality, we choose the LOS in question to 
be along x — 0. Also, for simplicity, we assume that there is no 
variation in B ± (Q, z) over the small portion of the line of sight in 
which we are interested, i.e. 



B ± (Q,z) ~ B ± (Q,Q) = const. 



(9) 



With the magnetic fields specified, we can now calculate F. 
From Eqs. ©, ©, ©, and ©, we find that 



F(4>) 



S 



dz 



axB'mz 2 



a s \B ± (0)\ 2 e 2i * 



(10) 



We recall that we assumed that B ± has a fixed orientation, which 
implies that^- is also fixed. 

To perform the z integration, we use the property of the delta 
function 

£ (id 

x:g(x)=Q 5 v ' 

where x are the roots of g{x), the sum is over all roots, and the 
prime indicates differentiation with respect to x. The argument 
of the delta function has two roots, z - + ^2(<p - (po)l(a\B' z {Q)). 
This means that the polarized emission from two values of z will 
contribute to F at a single value of <p. After performing the z 
integration, we find that 




Fig. 1. The normalized Faraday spectrum, F((p - (po)/D, for a 
caustic. The solid and dashed lines apply when B' 3 (0) > and 
Bj(0) < 0, respectively. 



F(0) = 



D 0[(0 - fa)B' 3 (0)] 



where 



D = 



V0 ~ 4>o 



2a s \B ± (0)\ 2 



(12) 



(13) 



and is the Heaviside function. 

The final form of F, as shown in Fig. Q~l is as we expect. 
The spectrum diverges at = (po, and is asymmetric with a tail 
extending to one side or the other of <po. The direction of the tail 
is determined by the sign of B^O), as indicated by the Heaviside 
function in Eq.fT2l 

Since in many cases caustics will be unresolved in obser- 
vations (see Sect. [3]), we should consider the integrated inten- 
sity over the range 6(p, which is taken to be approximately the 
FWHM of the point spread function in <p (often called the rota- 
tion measure spread function or RMSF). The integrated inten- 
sity, f, is 



r 



d<p F(cp) 



Together with Eqn.Q~2] this gives 



T = \2D V^l = 4a. s |B_l(0)| 2 



\flitf 3 (0)" 



(14) 



(15) 



As expected, we find that small values of B' 3 lead to strong caus- 
tic features. Large B ± values will also result in large values of T ■ 
However, since the B-field strength is unlikely to deviate signif- 
icantly from the average and the 1 / ^JW^ term can be arbitrarily 
large, we argue that strong caustics will mostly be created in 
regions where B' 3 is small. A result of this effect is that large 
"spikes" in the Faraday spectrum can result from intrinsically 
weak emission regions in physical space. 

We have only considered an illustrative case above, assum- 
ing constant thermal and relativistic electron fields, and the sim- 
plest possible magnetic field that includes a LOS reversal. More 
complicated fields will not have a significant effect on the fun- 
damental structure of Faraday caustics. Regardless of the form 
of the LOS magnetic field reversal, a caustic will appear as 
a strongly peaked, asymmetric source in Faraday space. The 
higher order terms in the expansion of B3, and the variation in the 
particle and B ± fields, will affect the specific shape of the "tail." 
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Equation [12] is the asymptotic form of F for a caustic when ap- 
proaching <pQ. The range of (p over which this expression applies 
depends on how closely the real fields (magnetic and density 
fields) are described by our assumed fields. 

The precise functional form of F{(f>) will also dictate the inte- 
grated flux over the feature. If F(<p) falls significantly faster than 
-0 ' 5 , this could reduce the integrated flux to the point that the 
caustic may not be detectable. 

There is one scenario in which the fundamental description 
of a Faraday caustic will deviate from the one given above: the 
LOS magnetic field reduces to zero but does not change sign. In 
this case, B'^(Q) = 0, Bj,(z) °c z 2 , and (f>(z) oz z 3 - The Faraday 
spectrum will still be strongly peaked at <po, but the spike will 
be symmetric in this case, with tails extending outward on both 
sides of <po. 

It is also worth noting how a caustic appears in /1-space. 
From Eq. ([TJ, we see that the polarized intensity as a function 
of A, P(A), is found by taking the Fourier transform of Eq. (TTZV 
Doing so, we find 

PU) = De- 2iA2 *> i ~ 1 ~ i) ^ . (16) 
2/1 

Although this seems to be divergent for A — > 0, we recall that 
D contains a s , which contains the spectral component of the 
Faraday spectra, s(A). This spectral component ensures that the 
total physical spectra is finite. 

3. Observing Faraday caustics 

Now that we have characterized Faraday caustics (to first or- 
der), we can discuss the expected appearance of these features 
in observations. Since they occur whenever the magnetic field 
is aligned in the plane of the sky, we expect Faraday caustics to 
be a regularly occurring feature in RM synthesis observations of 
diffuse polarized fields. Depending on the instrument used and 
the details of the experimental setup, the appearance of Faraday 
caustics may vary, and as we later demonstrate, may not exhibit 
the characteristic asymmetry described in the previous section. 
Nevertheless, Faraday caustics will be present in many observa- 
tions and so it is important to know what to expect in order to 
properly interpret Faraday spectra and learn as much as possible 
about the underlying magnetic fields. 

In lieu of observational da t a, we m ake use of the Hammurabi 
software of IWaelkens et al.l (120091) to provide examples of 
Faraday caustics in a realistic Faraday spectrum. This soft- 
ware simulates the Galactic synchrotron emission based on user- 
defined models of the magnetic field, thermal, and c osmic-ray 
electro n distributions. Hammurabi was first used by ISun et al.l 
(l2008h to test Galactic magnetic field models by simulating 
the resulting polarized foreground emission in the Galactic 
plane. Other examples of its use may be found in iJansson et al.l 
d2009l) who compared images of the Galactic synchrotron emis- 
sion from WMAP5 with simulated emission using several ex- 
isting Galactic magneti c field models from the literature, in 
lOppermann et all (1201 Ol) who tested the novel LITMUS pro- 
cedure that atte mpts to detect m agnetic helicity in simulated 
sky maps, or in iJaffe et all (1201 Oh who developed parametrized 
magnetic field models in the Galactic plane. We modified 
Hammurabi to calculate the Faraday spectrum; at each step of 
the integration along a given LOS, we compute the Faraday 
depth and the intrinsic polarized emission from within the cell. 
This is stored in a three-dimensional (3D) cube with user-defined 
Faraday depth resolution. 







(a) 






J\J ~~ 



-40 -20 20 40 

<t> (rad/m 2 ) 





Distance (kpc) 

Fig. 2. A Hammurabi simulated Faraday spectrum and the asso- 
ciated Faraday depth and field parameters as a function of phys- 
ical depth, (a) The Faraday spectrum along the LOS at Galactic 
coordinates (90°, -15°) in arbitrary units, (b) Faraday depth in 
rad irT 2 . (c) The product of the LOS magnetic field and thermal 
electron density, which is proportional to the Faraday depth in 
code internal units, (d) The perpendicular magnetic field compo- 
nent squared, which is proportional to the synchrotron emissivity 
again in code internal units. 



For the examples presented here, Hammurabi wa s run us- 
ing the same simulation parameters as lSun et al 1 (120081) for their 
ASS+RING magnetic field model including a stochastic compo- 
nent (as well as the cosmic-r ay electron distribution a nd the ther- 
mal electron distribution of ICordes & Laziol (120021) ). except at 
lower angular resolution, and with higher LOS resolution (more 
integration steps). 

We find that the simulated Faraday spectra typically consist 
of some weak diffuse emission spread over tens of rad rrT 2 as 
well as many narrow, asymmetric "spikes". In Fig. |2(a)| we show 
a portion of the Faraday spectrum along the LOS at Galactic 
latitude and longitude (90°, -15°) that provides a clear example 
of the kind of features discussed in the previous section. Along 
this LOS, we find a prominent caustic feature at <p = -39 rad 
irT 2 , as well as another small cluster of a few caustic features at 
0=19 rad m~ 2 . The Faraday depth, B and n e values are shown 
as a function of LOS distance in Fig. |2(b)| The product of n e 
and B3 is related to the Faraday depth, and the quantity B\ is 
proportional to the synchrotron emissivity. The caustic feature 
at -39 rad irT 2 is associated with the field reversal at z ~ 2.2 
kpc, where the LOS magnetic field is zero and the synchrotron 
emissivity peaks. One of the caustics in the cluster around 19 rad 
irT 2 is due to the LOS magnetic field zero crossing at z ~ 1.6 
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-40 -35 -30 -25 



4> (rad/m 2 ) 

Fig. 3. A portion of the simulated Faraday spectrum from Fig. 
|2(a)| This is compared to the generic form for the Faraday spec- 
trum of a caustic, F{<p) oc 1 / ^ (see Eq. (TTZl i). shown as a 
dashed line. 



kpc, while the others are due to magnetic field structures that are 
not shown. 

The major caustic shown in Figure |2(a)| at <p = -39 rad 
rrT 2 clearly exhibits the primary features described in the pre- 
vious section. It is strongly spiked, has the characteristic asym- 
metric shape, and is reasonably well fit by the generic F((p) oc 
1 / - 0o form as shown in Fig. [3] The tail extends towards 
higher values of </> as a result of the positive slope of the Bt,(z) 
distribution at z = 2.2 kpc. 

It should be possible to resolve particularly strong and iso- 
lated caustic features such as the one shown here, thereby al- 
lowing one to observe the characteristic asymmetry and possibly 
even to reconstruct the shape of the tail of the spike. The abil- 
ity to resolve Faraday caustics would not only provide a means 
of identification, but also information about the magnetic field 
strength and structure. For instance, the direction of the "tail" of 
the spike is dictated by the sign of B' 3 (zo). If B3 changes from 
positive to negative, then the Faraday caustic will be extended 
when <p < <f>o. In addition, if one were able to fit the tail to Eq. 
( TTZl ). it would provide a measure of D and therefore of the ratio 

Bi(O)/ Vsp). 

Wha t are the necessary condition s for resolving a Faraday 
caustic? iBrentiens & de Bruvnl (120051) give the resolution, Stp, 
and maximum observable scale in Faraday depth space for an 
RM synthesis experiment as 

2V3 

*~aF« (17) 

max. scale — , (18) 

A 2 . 

mm 

where AA 2 is the sampled range in A 2 and A, mn is the mini- 
mum sampled wavelength. To reconstruct the shape of the caus- 
tic, its width must be several times Sip. Owing to their narrow 
shape, high Faraday depth resolutions will be required to resolve 
Faraday caustics. To achieve long "baselines" in /t 2 -space, an ex- 
periment must extend to low frequencies. This implies that low 
frequency instruments such as LOFAR will be ideally suited to 
observing these types of features. 

We also require that the maximum scale, given by Eqn. QjU 
is larger than the resolution, Sip, in order to have a chance of 
reconstructing the caustics shape. From Eqs. (IPTl and ([18}, we 
find that this condition is satisfied if 



>~ 1.5. (19) 

Vmin 

We consider the high-band LOFAR receivers that operate be- 
tween 1 10 and 250 MHz. The FWHM of the main peak of the 
point spread function in Faraday depth space is 6<p ~ 0.5 rad/m 2 
and the maximum <p scale is roughly three times larger. 

As an example of what we can expect caustics to look like 
in real data, we can "observe" our simulated Faraday spectra by 
evaluating its Fourier transform in /i 2 -space, sampling at some 
set of frequencies, and then Fourier transforming back to <p- 
space. Figures |4(a)| and |4(b)| show the Faraday spectrum of Fig. 
|2(a)| as "observed" using two different frequency ranges. The 
frequency range used for Fig. |4(a)| is 110-250 MHz, which is 
similar to that for the LOFAR high band antennas. The fre- 
quency range used for Fig. |4(b)| is from 300-900 MHz, which 
is similar to the low frequency portion of the GMIMS survey 
(IWolleben et al. 2009). We include 512 frequency samples over 
each range. Most next generation instruments will be able to ob- 
serve using many more frequency channels than this. In each 
figure, we show the Faraday spectrum after deconvolution us- 
ing our own RMCLEAN s oftware that imple ments a procedure 
similar to that described bv lHeald et all (120091) . The RMCLEAN 
components are shown in red. 

In the "LOFAR" reconstructed spectrum, the stronger caus- 
tic is reasonably well recovered. The shape of the tail is not 
exactly as it is in the original spectrum, but the asymmetry is 
clearly visible. The small cluster of caustics at <p ~ 20 rad/m 2 
is resolved showing the three distinct peaks. The same cannot 
be said of the "GMIMS" spectrum. Owing to its lower resolu- 
tion (Sip ~ 5 rad/m 2 compared to ~ 0.5 rad/m 2 for the "LOFAR" 
observation), the individual peaks at </> ~ 20 rad/m 2 are not re- 
solved. The shape of the large caustic is completely lost, even 
though the tail of the feature extends over 15 rad/m 2 in the orig- 
inal spectrum. No asymmetry is present even in the distribution 
of CLEAN components. 

Observations with lower </>-space resolution, i.e. a smaller 
SA 2 , will of course fair worse than the above examples. We 
consider the Australian Sq uare Kilometre Arra y Pathfinder 
(ASKAP) POSSUM survey jGaensler et al]|2010l) that will op- 
erate over the range of frequencies from 700 MHz to 1.8 GHz. 
The 0-space resolution will be approximately 22 rad/m 2 , which 
is lower than the examples included here. While caustics will 
likely still be detected by the POSSUM survey, and other exper- 
iments operating at higher frequencies, they will simply appear 
as unstructured point sources. 

Having described the appearance of individual caustic fea- 
tures in an observed Faraday spectrum, we now investigate how 
these are distributed throughout the 3D Faraday spectral cube. 
Faraday caustics are generated at the boundaries of regions with 
opposing LOS magnetic field polarity. This boundary forms a 
surface in physical space. If the magnetic field distribution is 
well ordered, this surface will be smooth and map onto a con- 
nected surface in 0-space. Therefore, Faraday caustic surfaces 
will reflect boundaries between regions of opposite LOS mag- 
netic field polarity. 

As mentioned in the previous section, there are a symmetric 
variety of Faraday caustics that appear as a result of the LOS 
magnetic field approaching zero tangentially without changing 
sign. These features will appear when the boundary between re- 
gions of opposite LOS magnetic field polarity is aligned along 
the LOS, i.e. where a Faraday caustic surface folds over. 

In Figs. IB.ll and IB.21 shown in Appendix iBl we show six 
slices from a Hammurabi simulated Faraday spectrum at dif- 



6 



M.R. Bell et al.: Faraday caustics 



a 1 

8 
3 




(a) 


4- 
2 







i (rad/m 2 ) 



' 








(b) 


J 


1 h^l^X- 


Ihl J 


i. J lS? 





(rad/m 2 



Fig. 4. Mock observations of the Faraday spectrum shown in 
Fig. |2(a)| The spectrum has been observed over two frequency 
ranges: (a) 110 to 250 MHz, similar to the LOFAR high band 
antennas and (b) 300 to 900 MHz similar to the low frequency 
portion of the GMIMS survey. In each figure, the Faraday spec- 
trum has been deconvolved using RMCLEAN. The black line 
represents the restored spectrum, while the RMCLEAN compo- 
nents are shown as vertical grey lines (red in the online version). 



ferent Galactic latitudes. The top panel of each slice shows 
the Faraday spectrum as a function of Galactic longitude and 
Faraday depth, while the bottom panel shows the LOS magnetic 
field strength as a function of Galactic longitude and LOS dis- 
tance. The scale of the magnetic field strength has been com- 
pressed to highlight the regions of opposite polarity. White re- 
gions indicate a negative LOS magnetic field, black regions a 
positive field, and the colored sections highlight the boundaries 
where B3 « 0. 

In each slice, we can see that the smoothly connected bound- 
aries between regions of opposite LOS magnetic field polarity 
produce connected lines of Faraday caustics in 0-space. The 
caustics can be seen clearly in the Faraday spectrum, appear- 
ing as sharp, bright lines with faint tails extending to one side. 
Watching how the lines move from frame to frame, we can get 
a sense of how the caustics are distributed in 3D, appearing as a 
surface of bright polarized emission. In Figure IB. II at Galactic 
latitude -117.6°, a continuous line of caustic features extends 
from Galactic longitude 15° and <f> « -5 rad irT 2 to Galactic lon- 
gitude 60° and 0^0 rad irT 2 . This surface reflects the B3 = 
boundary at z « 0.2 kpc. This continuous line remains dominant 
until frame c, where the polarized emission has shifted to <p = 
rad irT 2 . At this point, a second B3 = surface, where the LOS 
magnetic field changes from negative to positive, is forming be- 
tween 4 and 6 kpc. The B3 = surface between longitudes -15° 
and 15° and at z ~ 0.6 kpc produces the line of Faraday caustics 
between <p ~ and 5 rad irT 2 . 



B 

z 




Fig. 5. A diagram illustrating how a small-scale fluctuation about 
a larger scale magnetic field variation may split a single Faraday 
caustic (occurring at the zero crossings) into several smaller 
caustics. An observation will not be sensitive to fluctuations be- 
low a given scale length that depends on the resolution of the 
experiment. 

In Figure IB. 21 the rightmost B3 = surface produces the 
dominant lines of Faraday caustics extending to positive <p values 
in the Faraday spectra. The Faraday caustic surface associated 
with the leftmost B3 = surface is now at <p — rad/m 2 . We 
note that the tails of the caustic spikes extend between the two 
surfaces since the LOS magnetic field changes from negative to 
positive causing one set of caustics, and from positive to negative 
for the other. 

Depending on the smoothness of the magnetic field varia- 
tion, and the resolution of a given experiment, the caustic sur- 
face may consist of a single, well-defined peak or a tight cluster 
of smaller spikes. We draw attention to the LOS magnetic field 
sketched in Fig. [5] where the black line represents a large-scale 
variation in the magnetic field, while the grey curve (red in the 
online version) represents small-scale fluctuations about this dis- 
tribution. For a smooth variation such as that represented by the 
black line, a single Faraday caustic will be present. In the case of 
the red curve, three caustics will be created close to one another 
in </>-space. The limited resolution of a given experiment acts to 
effectively smooth the magnetic field distribution by setting a 
minimum length scale below which the experiment is insensi- 
tive to variations. In other words, the §<p resolution also ensures 
a finite resolution in z-space. Hence, we expect Faraday caustic 
surfaces to be either sharp, smoothly connected features in the 
Faraday spectra, or a more fractured boundary depending on the 
scale of fluctuation of the magnetic field and the resolution, 5<p. 

4. Distribution of Faraday caustics in a Gaussian 
random magnetic field 

We have just shown how Faraday caustics can provide informa- 
tion about the structure of the LOS magnetic field from their 
shape and 3D distribution. It is variation and turbulence in the 
magnetic field that leads to caustics; after all, there would be no 
such features in a purely uniform magnetic field. It may be possi- 
ble to recover information about the statistical properties of the 
underlying magnetic field in a source by studying the distribu- 
tion of caustics within its Faraday spectrum. We therefore now 
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investigate the production of Faraday caustics by a Gaussian ran- 
dom magnetic field and calculate their probability distribution as 
a function of strength, i.e. the luminosity function. For detailed 
derivations of all results in this section, we refer to AppendixlAl 
We assume that the magnetic field distribution is a Gaussian 
random field with given covariance matrix M. For a Gaussian 
random field, the covariance matrix is also equal to the corre- 
lation tensor. The translationally invariant magnetic correlation 
tensor for homogeneous and isotropic magnetic turbulence is 



M^) = (B;(x)Bj(x + r)) 
Mij(r) = M N (r)6 u + (M L (r) - M N (r))^ 
+ M H (r)eij m r m in real space 

M^k) = M N (k)(6 U - y-) 

k 

- ieij m H(k)-^- in Fourier space (20) 



dSubramanianl 1 1 999l e.g.) where the longitudinal, normal, and 
helical auto-correlation functions, Mi(r), M^(r), and M#(r), re- 
spectively, only depend on the magnitude of the separation be- 
tween two points in space. We note that M^(k) can be expressed 
in terms of the one-dimensional (ID) magnetic field power spec- 
trum efl(fc) as 



e B (k) 



k 2 M N (k) 
8^ ' 



(21) 



The probability of measuring an integrated intensity of T at 
some location zo, given that there is a Faraday caustic at zo and 
given the magnetic covariance matrix M, is 



P(lF|caustic, M) 



P(T, caustic|M) 
P(caustic|M) 



(22) 



The denominator, P(caustic|M), is fairly straightforward to com- 
pute 

P(caustic|M) = I DB 6(B 3 (z)) 0{B, M), (23) 



where DB indicates a phase space integral over all possible B- 
field configurations. In addition, Q{B, M), a Gaussian probabil- 
ity density function over this phase space, is 



0(B,M) 



1 



(24) 



where vertical bars indicate that the determinant is to be taken. 
Here, the notation A J B implies a multidimensional scalar prod- 
uct such that 

The sum is over the field components and the integral is over 
physical space. 

We find that the probability of a caustic occurring (i.e. that 
Biizo) = 0) is 



P(caustic|M) = 



In 



M 33 (0)' 



(26) 



where M33 is the component of M that describes the two-point 
correlation between z components of the magnetic field. 



The calculation of the numerator is significantly more in- 
volved and as we describe in Appendix [A] Eq. (l22l i can only be 
computed analytically in the limit where T is large (we define 
this more clearly below). In this limit, we obtain the result 



P(;F|caustic, M) 



448a^0 3 



Ml(0) 



^27t|m;' 3 (0)| 



(27) 



We find that, over a limited range in f, the distribution of the 
integrated flux of caustics will be proportional to . We note 
that this indicates that the number of caustics per logarithmic 
interval of T is oc ^" _2 and thus that the Faraday caustic flux per 
logarithmic interval of T goes as T~ l , which is finite as T — » 00. 

The fore-factors in Eq. (l27l depend on the properties of the 
magnetic field and the turbulent flow of the fluid in which they 
are situated. The term Mn(Q) is simply the total magnetic energy 
density. We also introduce the Taylor microscale, It, which is 
defined by 



d 2 C(z) 



dz 2 



(28) 



(see e.g. iTennekes & Lumlevll 19721 sec. 6.4), where C(z) is the 
covariance in z normalized to unity at z = 0. This length scale 
represents the largest scale on which dissipation is important in 
a turbulent flow. The second derivative of the LOS component of 
the magnetic correlation tensor, M" 3 (0), is proportional to £b/I t - 
We can therefore rewrite Eq. (l27t as 



PCF|caustic, M) = 



448a? 



» 3 el% 



(29) 



where e# is the total magnetic energy density. 

By measuring the luminosity function of Faraday caus- 
tics in the regime, we can see that one is in principle 
able to measure the Taylor microscale in a turbulent medium. 
iFletcher & Shukurovl (120071) show that this length scale is also 
measurable by observing the mean separation between so-called 
"depolarization canals," which are lines of zero polarized in- 
tensity in a diffuse polarized field. These canals have been ob- 
served in maps of dif f use polarized Ga l actic emission by, e.g. 
lUvaniker et all (1 19981). l Haverkorn et al.l (l2000l) . iGaensler et all 
d2001h . and iReich et alj d2004l) . Observations of depolarization 
canals and Faraday caustics in tandem may be a powerful tool 
for studying the turbulent properties in, e.g. the ISM. 

A rather important caveat, applicable to both Faraday caus- 
tics and depolarization canals, is that without sufficiently high 
resolution one will not be able to measure this scale. In the case 
of caustics, as discussed in the previous section, the effect of 
small-scale fluctuations is to split a single Faraday caustic into 
a tight bundle of spikes. Without enough resolution, an observer 
may simply count such a bundle as one large spike. This has the 
effect of altering the normalization of the number distribution of 
caustics in f, which would lead to an incorrect measurement of 
It. In principle, the same effect applies to depolarization canals, 
where what appears to be a single canal of net zero polarized in- 
tensity may be a fine network of such canals. As described pre- 
viously, the effect of finite resolution is to smooth the quantity 
being measured, e.g. the LOS magnetic field distribution in the 
case of Faraday caustics, thereby throwing away information on 
any smaller scales. If the Taylor microscale is smaller than the 
smoothing scale set by our experimental resolution, then it will 
not be measurable. We suspect that this may be the reason why 
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the Taylor scale of the Milky Way is meas ured to be much larger 
than ex pected in the example presented bv lFletcher & Shukurovl 
d2007h . 

The result given in Eq. (l29l l is valid for strong spikes, i.e. 
when B' 3 (0) is small, but we expect that below some value of 
T = Tow the distribution will flatten. This is required to ensure 
that the integral of P(:F|caustic, M) is finite for T — > 0, but we 
can also see that this will be the case by considering the approx- 
imation applied during our derivation. 

To permit us to compute the probability analytically, 
we neglect an exponential term in Eq. ( |A.9t . This term is 
exp(-fiKB 4 L /'F), where K = \ba 2 64>lci\. This is, of course, only 
valid if the argument of the exponent is much less than one. The 
approximation breaks down when this term approaches unity, 
which leads to the condition for Ti ow 



where 



f 

Jo 



dke B (k) 



(30) 
(31) 

(32) 



is the average magnetic energy density, e B (k) is the ID magnetic 
power spectrum, and 



poo 

n J dke B (k)/k 
JT°° dke B (k) 



(33) 



is the magnetic correlation length. 

If we wish to evaluate Eq. jA.9i in all regimes, we can no 
longer neglect this term. While this full calculation is beyond the 
scope of this paper, we can consider what effects the inclusion of 
this term might have on the result. To proceed without approxi- 
mations, we would use a pertubative expansion in Feynman dia- 
grams as is used in quantum field theory. The first term in such 
an expansion would be negative and oc T~ 2 , the second positive 
and oc T -1 , and so on. The summation of these terms results in 
a turnover of the spectrum near Tiow 

We also expect the distribution to steepen at larger values of 
T because of depolarization. As discussed above, large values 
of T imply that the LOS magnetic field changes over long dis- 
tances. As this happens, it is increasingly more likely that the 
position angle of the polarized emission at the values of z that 
contribute to the polarized intensity at a single value of <p will be 
uncorrected resulting in a reduction in the integrated intensity 
because of depolarization. This reduces the number of caustic 
features at high values of T. 

To investigate the conditions for a steepening of the proba- 
bility distribution of caustics, we consider the correlation of the 
sky plane magnetic field as a function of LOS distance 



(b 2 A0)} 



(34) 



where the average is over all possible field configurations. The 
procedure for computing these expectation values is similar to 
that used in Appendix|A]to compute P{T\caustic, M) and there- 
fore not included here. The result is 



'R 



3M N (2z) 7M 2 N {z) 



IOMjvCO) 10Af2(0)' 



(35) 




Fig. 6. The normalized correlation of the sky plane magnetic 
field as a function of koz for a few values of the power law index 
a. Dashed red a = 4/3, solid black a = 5/3, dot-dashed blue 
a = 7/3. 



A small value of ft indicates that, on average, the sky plane mag- 
netic fields are uncorrected. When this is true, the luminosity 
function is steeper than T~ 3 ■ 

As an example of estimating Tiow an d ^> we model the mag- 
netic power spectrum as a broken power law 



e B (k) = e Q | — 



"5 



(36) 



For such a power spectrum, the luminosity function turns over at 



Tio„ = 7?<K 



eok Q ' 



a - I b+ 1 r la b\ '- 



2' 2, 



(37) 



where S(x) is the Beta function. 
We can also compute % because 



M N (z) 



I 



dkM N (k)e 



-ikz 



and 



M N (k) 



8n 3 e B (k) 
k 2 ' 



(38) 



(39) 



The resulting form of 'R{k( l z) is shown for a few representative 
values of a in Fig. [6] For a Kolmogorov spectrum, a = 5/3. In 
each case, b — 2. We note that % is independent of the parameter 
eo- 

A calculation of the precise shape of the distribution for low 
and high values of T is beyond the scope of this paper. A sketch 
summarizing the expected shape of the probability distribution, 
represented as the number of caustics per logarithmic scale of 
T as a function of log T, is shown in Fig. [7] The T~ 2 regime 
is valid given the assumptions outlined above for sufficiently 
strong caustics and until depolarization becomes important, as 
described above. We expect that the exact shape of the distribu- 
tion at and below Tiow will be more strongly dependent on the 
magnetic field statistics because weak caustic features depend 
more strongly on |Z? ± (0)|. 

In summary, we have calculated the probability density func- 
tion of caustics with a particular strength in a Gaussian random 
magnetic field. This result illustrates the type of analysis that is 
possible and represents the simplest reference case for future ob- 
servations of the luminosity function. We caution that the result 
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log T 



Fig. 7. A sketch of the number distribution of caustic spikes per 
logarithmic interval of T showing the limits to the T^ 1 regime 
and the rough shape of the distribution beyond these limits. At 
high values of T, the distribution steepens owing to depolariza- 
tion effects, while at lower T values the distribution turns over 
because the total number of caustics is finite. We refer to the text 
for more details. 

has a limited scope for several reasons. For one, we have as- 
sumed an unrealistic magnetic field distribution. This was done 
for simplicity, but also because if one is unaware of any statis- 
tical properties of the magnetic field other than the two point 
correlation, the assumption of a Gaussian field is the most ap- 
propriate starting point. Any further knowledge would always 
be included by expanding around the Gaussian case. While we 
expect our result in the regime to be only weakly depen- 
dent on the magnetic field statistics, realistic magnetic field dis- 
tributions are likely to be significantly non-Gaussian. We have 
not investigated the degree to which this will effect our results. 
Secondly, the result in Eq. ( |29] l is only valid over a limited range 
of T values. The precise extent of this range is unclear at the mo- 
ment, although we have presented some estimates. It is possible 
that the flattening and steepening regimes that we describe are 
near one another or even overlapping. A more detailed inves- 
tigation into the flattening regime, for example, would require 
much more complicated calculations and, because this regime 
depends strongly on the magnetic field statistics, a more realistic 
model. 

Lastly, Eq. d29b is truly only valid for very high resolution 
observations or simulated data. As discussed above, the result 
will change depending on 0-space resolution. In our treatment, 
we simply count the number of zero-point crossings of the LOS 
magnetic field. In observations, two or more of these events in 
physical space may occur at nearly the same (f> location and may 
therefore be counted as a single caustic. This would have the ef- 
fect of changing the overall normalization of the distribution (the 
total number of caustics would be reduced), but the predic- 
tion should still be valid. Some additional analysis of simulations 
would be helpful in assessing the extent of this effect. 

5. Discussion and conclusions 

We have introduced sharp, asymmetric features in the Faraday 
spectrum of a diffuse polarized source that appear as a result of 
a change in the direction of the LOS magnetic field. We have 
called these features Faraday caustics, and derived their funda- 
mental properties by considering a simple magnetic field distri- 



bution near the point where B3 = 0. The general properties of a 
Faraday caustic are: 

- They are singularities in the Faraday spectrum. The polarized 
intensity peaks sharply at the Faraday depth 0o that maps to 
the physical depth zo, where £3(20) = 0. 

- The Faraday spectrum extends only to one side of the sin- 
gularity. The direction in which the tail extends depends on 
the slope of the B 3 (z) distribution. If B' 3 (z) > 0, then the 
tail extends to <p > 0o- An exception occurs when Z?3(z) 
approaches zero tangentially without changing sign. In this 
case, the caustic is still strongly peaked, but symmetric. 

- Strong caustics will predominantly be created when the LOS 
magnetic field changes over large physical distances i.e. 
when BUzo) is small. In this case, the polarized emission 
from even a relatively weakly emitting region (in physical 
space) can appear as bright point sources in the Faraday 
spectrum because the emission over a long LOS distance can 
pile up over a narrow range of Faraday depths. 

- In a 3D Faraday spectrum, Faraday caustics will appear as 
continuous sheets or surfaces. These surfaces are related to 
boundaries between regions of opposite LOS magnetic field 
polarity. 

- The effect of small-scale fluctuations about larger-scale vari- 
ations will be to split a single large Faraday caustic into a se- 
ries of smaller, closely packed caustics. Without sufficiently 
high resolution, this splitting will not be observed. In this 
way, the finite resolution of an experiment acts to smooth 
the Bt,(z) distribution. 

We find that important information about the magnetic 
field distribution can be recovered by observing these features. 
Reversals of the LOS magnetic field can be identified and if the 
asymmetry of the Faraday caustics is recovered, one will even be 
able to indicate in which direction the reversal occurs. We have 
also shown that there may be great value in studying the statisti- 
cal properties of Faraday caustics. We have computed the prob- 
ability distribution of the polarized intensity of Faraday caustics 
in a Gaussian random magnetic field. The number of caustics 
per unit T , which is the integrated intensity of a single caustic, 
goes as . There is a lower T limit at which the distribution 
flattens and turns over, and a high T limit where the distribution 
steepens because of depolarization effects. More work is needed 
to more clearly define these limits, but this should be a fruit- 
ful endeavor because we expect that these limits will depend on 
the statistical properties of the magnetic field distribution. We 
show that if one has sufficiently high resolution and measures 
this distribution in the regime one can recover the Taylor 
microscale. This scale is the largest distance over which dissipa- 
tion is important in a turbulent flow. 

Faraday caustics will not be the only source of point-like 
structure in the Faraday spectrum of a diffuse source. Point 
sources in Faraday space will also occur as a result of a dis- 
crete source along the LOS, such as a pulsar or any strong, com- 
pact emitter. These discrete sources will also appear as compact 
sources in the sky plane and will often have associated structure 
in total intensity making them easy to distinguish from Faraday 
caustics. If a caustic is resolved, the asymmetric shape will of 
course be an obvious means of identification. In the case of struc- 
tures more extended in Faraday depth space than the largest mea- 
sured scale (c.f. Eq. (TT8ll). the extended nature of the source will 
not be recovered and the ed ges of the structure will appear to be 
point-like (see e.g. Fig. 2 in lFricket ail (1201 Pi) - ) . 

It should be noted that caustics may already have been ob- 
served. Though we have not analyzed the data directly, the 
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publi s hed data from the Perseus field (Ide Bruvn & Brentjensl 
l2005t lBrentjensll2010t) contains many point-like features, some 
of which seem to be connected in what may be the Faraday caus- 
tic sheets we discuss above. A re-analysis of these data may pro- 
vide the first observational evidence of Faraday caustics. 

New observational techniques and imaging algorithms may 
be needed in order to optimally observe Faraday caustics. When 
investigating t he effectiveness of a wavelet-based approach to 
RM synthesis, iFrick et ail (1201 ll) noted sharp, point-like fea- 
tures in a Faraday spectrum produced by the small-scale fluctu- 
ations included in their magnetic field distribution model. They 
also first showed the characteristic asymmetric caustic profile 
resulting from a large-scale magnetic field reversal in their Fig. 
5. Neither of these caustic features were discussed in detail, or 
identified as caustics, but the authors showed that wavelet-based 
RM synthesis is effective at separating small and large scale fea- 
tures in the Faraday s pectrum and as a re sult the wavelet-based 
algorithm outlined in IFrick et all (1201 Ol) may be a useful tool 
for observing Faraday caustics. In addition, matched filter algo- 
rithms could be helpful if one wants to more accurately recon- 
struct the shape of the "tails", which contain most of the total 
polarized flux and information about the magnetic field distribu- 
tion. 
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Appendix A: Derivation of the probability distribution for caustics 

We derive the probability of measuring a spike with an integrated intensity f given that the spike is due to a caustic (i.e. B 3 = 0) 
and a magnetic covariance matrix M. We limit our discussion to the regime where Eq. (ITSt adequately describes the integrated flux 
of the caustic. 

In Sect. |2j we found that the integrated polarized intensity of a caustic is T — 4a s -\[6^ \B ± (0)\ 2 / -^B'^tyai. We assume a 
Gaussian magnetic field distribution, Q{B, M). The probability distribution of integrated intensities T, given that Bj = and the 
magnetic correlation tensor is M is 

_ P(T, causticlM) 

PCFIcaustic, M) = — — — . (A.l) 

P(caustic|M) 

The denominator, which is the probability that a caustic occurs in a Gaussian random magnetic field, is relatively straightforward to 
compute. We impose the condition that B 3 (Q) = using a delta function in order to compute the probability 

P(caustic|M) = j DB S(B 3 (0)) Q(B, M), (A.2) 



where the integral is a path-integral over all possible realizations of the magnetic field. We then replace the delta function using the 
Fourier transform 

P(caustic|M) = -I f DB f dr] e 2 *^ 3 e -W*- 1 *. (A.3) 

V|27rM| J J 

We now replace r\B 3 with rjda B and after completing the square we have 



P(caustic|M) =^-^j j J drj exp - ^ {B - J rt Af) t M~ l {B - /ikf) + ^MJ 



(A.4) 



where /' = 2mri8i Z 8{z). The B integral is the integral of a Gaussian function, which gives a factor of y/\2nM\. The r\ integral is also 
a Gaussian integral. After integration, the result is 



P(caustic|M) = 



2n 



(A.5) 



The calculation of the numerator of Eq. (IA.U is more complicated. We again start by imposing our conditions using delta 
functions 



P(T, caustic|M) = \ DB 6(B 3 (0))6(r 



a s ^B 2 ± (0) 



^f\2lri\4\ 

Using Eq. ( TTTb . we rewrite the second delta function in terms of B 3 



1 



V# 3 (0)ai ; 

/ 4a,B 2 {0)yM\ B 4 I , B 4 .(0)\ 
. luy y = 4K ^L d l B > (o) _ K-^ , 

v -jB'(0)ai I r 3 \ 3 r 2 ) 



B^M l B 



where K = \6a 2 6(/>/a l . 
Inserting this into Eq. 



we now proceed with the calculation using the Fourier representation of the delta function 



(A.6) 



(A.7) 



DB 6d> 2 B 4 
6(B 3 (0))2K 



fs(B' 3 (0)-K-±>) 



Vi2lrMi 
2K6(p 2 
a 2 T 3 



a\T 3 



1 .,. , 
-B ] M~ l B 



J J 2nJ 2n ±K ' F 



exp[/77Z?3(0)] exp 



/ , KB 4 (0)\ 



(A.8) 



We introduce two generalized fields yi = 8a 6{x - xq) 6(y - yo) /a and 77' = <5o 6(x - xo) S(y - yo) r\ that permit us to work with 
the total magnetic field B(z) rather than the LOS-component B 3 (z) in Eq. (IA.8b : 



4K r DB C da Cdri , en% \ /«"5*(0)\ 



exp 



- -B^M~ l B + i(rf + ytdjfB 



(A.9) 



The functional integral over DB cannot be solved analytically because of the B^(z)-term in the exponential. In general, we need to 
employ diagrammatic perturbation theory in order to proceed. 

If we restrict our attention to strong caustics, we can proceed analytically. This is a reasonable restriction because it will be the 
strong caustics that are the most likely features to be observed. We note that a strong caustic can occur when either B ± is large 
or B' 3 is small. Strong caustics caused by an exceptionally large sky plane magnetic field are unlikely because of our assumption 
of Gaussian statistics, and therefore we concentrate on caustics with large T owing to a slowly changing LOS magnetic field. For 
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these caustics, we can neglect the B^(z)/^ r4 -term and solve the remaining problem by introducing a generating functional J and 
repeatedly applying the Gaussian integrations 



4K f DB r du r dri , , , ^ 

—7 , ^- ^ (B?(0) + B\ (0))- exp 

"F 3 J ^OT] J 2^J 2w|U/i(z)J + 

(_L_f + (_L_) 2 f f^exp 
U/ifcV \6J 2 (Z)> \ J 2jtJ 2n F 

/ g \ 2 I 6 \ 2 f fdp_ Cdv 
Uy,(z)/ + U/ 2 (z)/J J 2tt J 2^ CXP 

+ ( - \(i + I?M(p!d 3 ) - -Qi'd 3 ^M(r]' + J)) +17]'* M J 



- i^AT'B + i(7j' + fi'd^B 



4K 

F 3 

4K 

F 3 



-B f M~ l B + (vrj' + ifi'd 3 + jfB 



-{ir}' + ifi'di + JyM(iT]' + ifi'd 3 + J) 



- \rpMrf - i^'5 3 ) 1 "M( yU '5 3 ) + l -fMJ 



,/=o 



,/=o 



=0,because M^{z,z') is at its maximum for z,z f =Q 



4K 



1 



y/2n\M£(0)\ 
4K 1 



./=o 



=/ 



!F 3 V27r|M 33 (0)| V2^|M^(0)| 



4# 



1 



1 



5" 3 V2tt|M 33 (0)| V^WSM 
28a: 1 1 



f3(M 2 !(0) + M 2 2 (0)) + M 22 (0)M„(0) + (M 4 l3 (0) + M 4 23 (0) + 2M 2 3 (0)M 2 3 (0))J 



1 t 1 i l 

-J ] MJ 

2 2 M 33 (0) 



'33 



Mi, = 



28A: 



1 



1 



:F 3 V2tt|M 33 (0)| V 2 ^l^ 33 (0)1" 
2SK 1 1 



T 3 V2tt|M 33 (0)| V2^' 3 (0)| 



!F 3 V27r|M 33 (0)| V2^|M 3 ' 3 (0)| 

2 



r </& 3 
J (&o 



M N (k) 



dk 



3 £B(k) 



k 2 



(A.10) 



In the last step, we identified the covariance matrix M with the magnetic correlation tensor given by Eqs. d20l and (|2TT i. We used 
the property Afy(0) = Mai(0) <5,j. Finally, from this result and Eq. (|26| |. we conclude that in the bright caustic limit 



P CF|caustic, M) 



P (T, caustic|M) 448a 2 ^ 2 M 2 (0) 



P(caustic|M) r 3 ai ^2tt\M^U)\ 



448a 2 (50 3 



T 3 a\ y/27T\M' 3 ' 3 (Q)\ 



dk 



e B (k) 



k 2 



2 = 448a 2 r50 3 [/^ffif 

r 3 a 2 ^y dk 3 €B{k) ^ : 



(AM) 
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Fig. B.l. The first three slices of a Faraday spectrum cube of a Hammurabi simulation and the associated LOS magnetic field 
strength. Each pair shows a frame at a different Galactic latitude. Top: A grey scale image of polarized intensity (in arbitrary units) 
as a function of Faraday depth and Galactic longitude. Bottom: The LOS magnetic field strength as a function of LOS distance 
(in pc) and Galactic longitude. The scale has been compressed to highlight regions with opposite polarity. Black indicates positive 
magnetic field strengths and white negative. The colored area is where B3 as (appears as grey scale in text version). 
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Fig. B.2. Three more slices of a Faraday spectrum cube of a Hammurabi simulation, following those in Fig. IB. II 



